Quantum spin correlations in an organometallic alternating sign chain 
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High resolution inelastic neutron scattering is used to study excitations in the organometallic mag- 
net DMACuCls. The correct magnetic Hamiltonian describing this material has been debated for 
many years. Combined with high field bulk magnetization and susceptibility studies, the new results 
imply that DMACuCls is a realization of the S = 1/2 alternating antiferromagnetic-ferromagnetic 
(AFM-FM) chain. Coupled-cluster calculations are used to derive exchange parameters, showing 
that the AFM and FM interactions have nearly the same strength. Analysis of the scattering in- 
tensities shows clear evidence for inter-dimer spin correlations, in contrast to existing results for 
conventional alternating chains. The results are discussed in the context of recent ideas concerning 
quantum entanglement. 

PACS numbers: 75.10.Jm, 75.40.Gb, 75.30.Et 



Low dimensional, dimerized quantum spin systems ex- 
hibit diverse physics ranging from simple spin gaps to 
exotic spin hquids[l|, 0, |3, A prototypical example is 
the 5 = 1/2 Heisenberg alternating chain (HAG), defined 
by the Hamiltonian = J2n Ji^'^ri-i-S2n + J2S2n-S2n+i- 
The HAG has been widely studied theoretically [^, @, 01 
and experimentally [1, Gonsidcring Ji > an- 

tiferromagnetic (AFM), ground state behavior of the 
HAG depends on the ratio a = J2/Ji- a = 1 corre- 
sponds to the quantum critical pure AFM chain. When 
I a I << 1 the system behaves as nearly independent 
dimers, and perturbation theory in a provides a straight- 
forward approach to calculate physical properties of the 
spin gap system. For a < 0, one has a realization of the 
antiferromagnetic- ferromagnetic (AFM-FM) HAG. For 
a « -1, the HAG maps onto the 5 = 1 AFM Hal- 
dane chain with the dimer viewed as a composite object 



lot I 111 . |12| |. The regime a « —1 is particularly inter- 



esting since the system exhibits intermediate behavior 
between dimer and Haldane physics Galculating 
physical properties must be done carefully since pertur- 
bation theory from the dimer limit is not a good approx- 
imation here. 

We report new inelastic neutron scattering (INS) re- 
sults for (GH3)2NH2GuGl3 [DMAGuGlg or MGGLl. We 
use state of the art linked-cluster calculations [H, [3| 
and bulk measurements to show DMAGuGls contains a 
realization of the AFM-FM HAG with a = -0.92(4). 
Moreover, in contrast to results on existing AFM-AFM 
HAGs [1, Q , there is evidence for spatially extended en- 
tanglement of the principle dimers via the FM coupled 
spins. 

DMAGuGla is monoclinic at room temperature with 



(a) AFM-FM alternating chain (b) AFM-AFM spin-ladder 
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FIG. 1: Cu^+ S = 1/2 sites in an ab plane of 
DMACuCls. Open and closed circles illustrate two distinct 
dimer bonds/chains in the LT phase [IE], (a) AFM-FM HAC 
scenario with Ji and J2 shown as solid lines and weaker ex- 
change, J', shown as dotted lines, (b) AFM-AFM SL scenario 
with J± and Jy . Lattice vectors correspond to HT monoclinic 
parameters [151 ]. 



(3 = 97.5° with copper-halide planes separated along c by 
methyl groups, such that magnetic coupling is only ex- 
pected in the a6-plane as shown in Fig.[ljl5,]. Thermody- 
namic and structural considerations were used to identify 
DMAGuCla as an AFM-FM HAG along the a-axis with 



1 [16|, llTl Il8l | , attracting interest since other ex- 



perimental realizations of AFM-FM HAGs have a <C — 1 
[l9l . [20j . Further measurements were interpreted as inde- 
pendent FM and AFM dimers with evidence for a tran- 
sition to long range order below T — 0.9 K 21, 2^- Pre- 



liminary INS in the disordered phase eliminated both in- 
dependent dimer and a-axis chain models [23] . Recently, 
it was shown that the low-temperature (LT), T < 285 K, 
structure is triclinic with two independent chains along 
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FIG. 2: T = 1.8 K constant Q scans of DMACuCla. (c) 
includes data at T = 40 K. Lines are based upon a global fit 
to the AFM-FM HAG model described in tiie text. Dashed 
lines are fitted backgrounds. 



the high-temperature (HT) &-axis [M |25|, cf. Fig. [T] 
Bond lengths in the LT structure favor HAC Cu-halide- 
halide-Cu or Cu-halide-Cu exchange paths, found to be 
relevant in other low-d magnets*2, J^, 4|. 

Deuterated samples were obtained by slow evaporation 
of a D2O solution of CuCla and (CD3)2ND-HC1. A 2.23 g 
single crystal mounted in the (hkO) plane was examined 
using the SPINS spectrometer at NIST configured with 
80' collimation before and after the sample and a cooled 
Be filter in the scattered beam. A flat PG(002) analyzer 
selected scattered neutrons at 5 meV. Constant wave- 
vector, Q, scans were performed at T = 1.8 K indexing Q 
in the HT monoclinic notation with LT lattice constants, 
a=12.05 A and 6=8.43 A^. 

Representative scans are shown in Fig. [2l Energy 
transfers were scanned up to hu; « 4.5 meV. A single 
peak disperses along k between 0.95(2) < Tilo < 1.66(2) 
meV as shown in Fig. [2ja)-(d). The mode is only weakly 
dispersive along h at the bottom (Fig.[2fe) and (f)) and 
top (Fig. [2Kg) and (h)) of the band with ss 0.2(1) meV 
dispersion along {h 1.5 0). Intensity reduction for tem- 
peratures large compared to the energy scales of exci- 
tations indicates a magnetic origin, c/. Fig. [UJc). The 
intrinsic peak width is small compared to the instrumen- 
tal resolution of « 0.25 meV FWHM. We do not see any 
reasonable evidence for gapless excitations or magnetic 
continuum scattering. 

We fit the constant Q scans with Gaussian peaks to 




FIG. 3: (a) Scattering intensity vs. hu along (OfcO). Solid 
lines are fits to dispersion based upon linked-cluster models. 
Dotted line is dispersion based upon perturbation theory for 
weakly coupled dimers. (b) Integrated intensity vs. Q. (c) 
Dispersion based upon Gaussian peak approximation. Lines 
in (b) and (c) are based upon a global fit to the AFM-FM 
HAG model described in the text. 



model (method) 


Ji or J± (meV) 


J2 or J|| (meV) 


AFM-FM HAG (E) 


1.406(8) 


-1.30(5) 


AFM-FM HAG (PT) 


1.319(8) 


-0.68(2) 


AFM-AFM SL (E) 


1.194(6) 


0.376(5) 


AFM-AFM SL (PT) 


1.319(8) 


0.34(1) 



TABLE I: Resulting exchange parameters from fitting dis- 
persion to AFM-FM HAG and AFM-AFM SL models using 
linked-cluster expansions (E) and perturbation theory (PT). 



determine dispersion and intensity variation in the (hkO) 
plane as summarized in Fig. [3l Lack of dispersion along 
h excludes the proposed a- axis as the chain axis, and 
significant dispersion along k excludes non-interacting 
dimer models. The k dispersion has periodicity 27r in 
reduced units, a quantum effect showing a lack of AFM 
order. The FM-FM HAC is ruled out by the absence of a 
second dispersive mode ]2^. Notably, the band max- 
imum [minimum] is at Q = 2rmT[Q = (2m + 1)tt] for 
integer m. This is unusual and rules out a HAC with 
a > 0. Models that can produce such dispersion include 
the AFM-FM HAC and the AFM-AFM spin-ladder (SL) 
with the HAC favored due to structural arguments. 
The dispersion for the HAC and SL models can be 
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calculated using perturbation theory, with an identical 
form to first order, nuj{Cl) = £1+62 cos(Q-u) (5l.l7l.li0l.l27j. 
where ei = Ji[J±] and £2 = ^-^[J\\] for the HAC [SL] 
models. It is more accurate to use state of the art linked- 
cluster expansion methods [l^, • Parameters obtained 
from fits of the data to these models are summarized in 
Tab.lll Dispersion alone can not differentiate between the 
HAC and SL models. Solid lines in Fig.[3^a) show the fit 
to the linked cluster calculation and the dotted line a fit 
to the perturbation theory. Perturbation theory is clearly 
inadequate to describe the AFM-FM HAC as evidenced 
by the large difference in fitting parameters relative to 
the linked cluster method. Although fits to the AFM- 
FM HAC dispersion characterize the magnitude of FM 
and AFM exchange, this comparison does not include 
structural information determining which bond in Fig. 
1(a) is FM or AFM. 

In principle, high field INS can be used to determine 
absolute values of exchange parameters assuming they 
are field independent [2^ . Unfortunately, to characterize 
DMACuCls accurately in this manner would require INS 
at magnetic fields well over hqH = 20 T; not possible at 
the present time. However, it is practical to examine high 
field bulk thermodynamic properties. Figure Ufa) shows 
pulsed field magnetization of DMACuCla up to 20 T at 
T = 1.6 K. The upper critical field where all spins are FM 
aligned is in the vicinity of 16 T. We also plot Quantum 
Monte-Carlo (QMC) calculations of the magnetization 
for 100 spins using exchange constants determined from 
cluster-expansion fits of hLu{Q) [2^. The AFM-FM HAC 
model reproduces the upper critical field indicating that 
the determined energy scales are more appropriate, but 
neither model accounts for the linear magnetization at 
low- field. 

Recall that the LT structure has two types of bond- 
ing, i.e. two sets of exchange parameters for quasi-ld 
chains along the monoclinic 6-axis. However, our INS 
measurements indicate only a single mode. The energy 
scale of the second set of quasi-ld chains is either much 
higher, lower or identical to that of the excitation we ob- 
serve. Higher energy excitations would manifest a higher 
energy gap in the spectrum, not seen in thermodynamic 
measurements. It is most reasonable that these excita- 
tions are at lower energies. The present experiment sets 
an upper limit of J' « 0.2 meV for this energy scale, in- 
distinguishable from the incoherent elastic scattering. In 
the absence of a measured exchange for the weakly cou- 
pled chains, we model these moments as free-spins. This 
is reasonable given the temperature scale of the magne- 
tization measurements. A conjecture of a free-spin con- 
tribution was also made in Ref. [s^l to account for mag- 
netization measurements. We therefore calculate M{H) 
as a composite model with half the moments included as 
free-spins. This accounts for the initial linear magneti- 
zation. As seen in figure 4(a), the composite AFM-FM 
HAC model calculation is a much better description of 



^oH (T) 
10 




150 200 
T(K) 

FIG. 4: (a) Pulsed field magnetization of DMACuCls powder 
at T = 1.6 K. Lines are QMC calculations based upon AFM- 
FM HAC and AFM-AFM SL models. Solid lines, labeled 
as composite chain or composite ladder, include a free-spin 
contribution. Symbols for data are plotted every 1000 points, 
(b) Single crystal magnetic susceptibility, x{T), and x{T)T 
for noH — 0.05 Tesla with H || a. Lines are the composite 
model AFM-FM HAC polynomial fit discussed in the text. 
Broken lines are the xiT) contribution of the AFM-FM HAC 
(dashed black) and free-spin (dash-dot green) components of 
the composite model, (c) X~^{T) based upon data in (b). 
Solid line is extrapolated linear fit described in the text. 



the the measurement than the composite AFM-AFM SL 
model. When M{H) is calculated using parameters de- 
rived from perturbation theory the situation is similar, 
but as expected the overall agreement with the data is 
not as good for either model. 

The calculations illustrated in Fig. Ufa) are performed 
with no adjustable parameters. The AFM-FM HAC 
model accounts very well for both the low field and high 
field limits with deviations in the vicinity of 5 to 10 T. 
This may be a consequence of a field-induced collective 
ordered phase i22] that is not accounted for in the calcu- 
lation. Willet et al. used finite chain calculations to ob- 
tain parameters that reproduce M[H) within a two chain 
model[2^. However, the exchange parameters they pro- 
pose are incompatible with the INS data reported here, 
and can be ruled out by the observation of only a single 
excitation. 

Figure lUJb) shows temperature dependent magnetic 
susceptibility, x{T)- The solid line is a fit using the 
AFM-FM HAC:free spin composite model, with x{T) for 
the HAC calculated using a polynomial representation 
derived from finite-sized scaling ^23|. The parameters 
Ji = 0.973(4), J2 = -1.23(5) meV and g = 2.096(1) are 
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reasonably consistent with INS results although clearly 
the latter provides a more direct measure of the exchange 
parameters. Figure [Ifc) shows x~^(T). A linear fit for 
120 < T < 270 K leads to a value for the Curie- Weiss 
temperature, 6 = —0.9{1)K. Linear fits to the compos- 
ite model calculations of x^^{T) yield 6 = — 0.4(3)if for 
the AFM-FM HAC, and 9 = -2.4(7)fs: for the AFM- 
AFM SL confirming the chain model is in better accord 
with the data. 

Neutron scattering intensities for a straight HAC can 
be calculated using the coupled-cluster expansions [l^ . 
Comparison to these using our HAC exchange param- 
eters and noting the broad wave-vector dependence of 
the integrated intensity, cf. Fig. [S^a) and (b), the short 
bond in Fig. [T] is likely AFM. However, the bonding ge- 
ometry in DMACuCla is complicated by in-plane zig-zag 
and small out of plane components. We incorporate bond 
vectors within a single mode approximation (SMA) 
to calculate the dynamic spin correlation function. The 
SMA is justified by the observation that the INS data 
show a single peak with no appreciable additional inten- 
sity, and that the calculated spectral weight is dominated 
by single particle excitations for a « —1 [12] . The SMA 
to the scattering intensity is 

%n{Cl,hij) cx ^M^^ Jd(So-Sd)[l-cos(Q.d)] 

d{hLu -huj{Q)). (1) 

(So • Sd) is the spin correlation function, d is a bond vec- 
tor, and the sum is over the AFM and FM bond, di and 
d2. We use 74 constant Q scans {huj > 0.5 meV) in a 
global fit to Eq.[T] convolved with the instrumental resolu- 
tion function and a dispersion based upon the AFM-FM 
HAC cluster-expansion using parameters in Tab. HI Pa- 
rameters include an overall scaling factor and the ratio 
A = (So • Sd2)/(So • Sdi). The results reproduce well 
the dispersion and intensity modulation and are shown 
as solid lines in Figs. [Hand [2Kb,c) with A = 0.26(4). Al- 
lowing exchange constants to be fit parameters results in 
Ji = 1.399(6) and J2 = -1.07(3) meV, and A = 0.27(4). 

The SMA has been applied to the AFM- AFM HAC, 
with the result that (So • Sdi) is consistent with the iso- 
lated dimer expectation of —3/4, and the empirical value 
for (So • Sda) = 00. The significance of this is dis- 
cussed by Brukner et al. f3l'|, who shows that the intra- 
dimer correlation is a manifestation of quantum entan- 
glement in the bulk system. Entanglement is localized 
to within one dimer for the AFM- AFM HAC. In con- 
trast for DMACuCls, the inter-dimer spin correlation is 
clearly non-zero. This is expected since as alpha tends 
to —00 the FM coupled spins can evolve towards a com- 
posite S = 1 entity, with a T = spin-spin correlation in 
the Haldane chain extending over several lattice spacings 
[3^ . Since the quantum correlation becomes spatially 
extended, adjacent dimers in the AFM-FM HAC can be 



thought of as more entangled than those in the AFM- 
AFM HAC. Notably, fitting the SMA for a purely AFM 
SL to the present data yields A « 0, indicating suppres- 
sion of interdimer entanglement for purely AFM coupled 
dimers in that case as well. 

In summary, high resolution INS, combined with 
high field magnetization and susceptibility measurements 
show that DMACuCls is a quasi-ld AFM-FM HAC with 
a « — 1, intermediate between weakly coupled dimer and 
Haldane regimes. The exchange parameters are deter- 
mined from the dispersion using the results of coupled 
cluster series expansions. Applying the SMA to analyze 
the scattering intensity shows that the inter-dimer spin 
correlation is significant, indicating that the AFM-FM 
HAC may be a model system for studying the effect of 
spatially extended quantum entanglement on bulk prop- 
erties. 
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